Multi-GeV neutrinos due to nn oscillation in Gamma-Ray Burst Fireballs 
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The long and short gamma-ray bursts are beheved to be produced due to collapse of massive 
stars and merger of compact binaries respectively. All these objects are rich in neutron and the jet 
outflow from these objects must have a neutron component in it. By postulating the nn oscillation 
in the gamma-ray burst fireball, we show that, 19-38 GeV neutrinos and anti- neutrinos can be 
produced due to annihilation of anti-neutrons with the background neutrons. These neutrinos and 
anti- neutrinos will be produced before the 5-10 GeV neutrinos due to dynamical decoupling of 
neutrons from the rest of the fireball. Observation of these neutrinos will shed more light on the 
nature of the GRB progenitors and also be a unique signature of physics beyond the standard model. 
A possible way of detecting these neutrinos in future is also discussed. 



I. INTRODUCTION 



Gamma Ray Bursts (GRBs) are flashes of non-thermal 
bursts of low energy (~ 100 KeV-1 MeV) photons and 
release about lO^-'^-lO^^ erg in a few seconds making them 
the most luminous object in the universe. They are 
known to occur at cosmo logical distance [11, 0. They fall 
into two classes; short-hard bursts (< 2 s) and long- 
soft bursts. The short bursts comprises about 25% of 
the total events and the rest are long-soft bursts. It is 
now widely accepted that long duration bursts are pro- 
duced due to the core collapse of massive stars the so 
called hypernovae. The origin of short-duration bursts 
are still a mystery, but recently there has been tremen- 
dous progress due to accurate localization of many short 
bursts by the Swift^, ^ and HETE-2[i] satellites. The 
afterglow observation of GRB 050709 at a redshift of 
0.1606[i| by HETE-2 and the Swift observation of af- 
terglow from GRB050709b at a redshift of 0.225[3j and 
GRB 050724 at a redshift of 0.2580 seems to support 
the coalescing of compact binaries as the progenitor for 
the short-hard bursts although definite conclusions can 
not be drawn at this stage. 

Irrespective of the nature of the progenitor, it is be- 
lieved that, gamma-ray emission arrises from the collision 
of different internal shocks (shells) due to relativistic out- 
flow from the source. A class of models call fireball model 
seems to explain the temporal structure of the bursts and 
the non-thermal nature of their spectra[ll, Sudden 
release of copious amount of 7 rays into a compact re- 
gion with a size cSt ~ 100-1000 kmQ creates an opaque 
7 — fireball due to the process 7 -I- 7-^ e"*" + e^. The 
average optical depth of this processQ is r-y-y ~ 10^^. 
This optical depth is very large and even if there are no 
pairs to begin with, they will form very rapidly and will 
Compton scatter lower energy photons. Because of the 
huge optical depth, photons can not escape freely. In the 
fireball the 7 and pairs will thermalize with a temper- 
ature of about 3-10 MeV. The fireball expands relativis- 
tically with a Lorentz factor F ~ 100 — 1000 under its 
own pressure and cools adiabatically due to the expan- 
sion. The radiation emerges freely to the inter galactic 



medium (ISM), when the optical depth is Tj~^ ~ 1. 

In addition to 7, pairs, fireballs may also contain 
some baryons, both from the progenitor and the sur- 
rounding medium. These baryons can be either free or 
in the form of nuclei. If the fireball temperature is high 
enough (more than 0.7 MeV), then it will be mostly in 
the form of neutrons and protons. The electrons associ- 
ated with the matter (baryons) can increase the opacity, 
hence delaying the process of emission of radiation and 
the baryons can be accelerated along with the fireball and 
convert part of the radiation energy into bulk kinetic en- 
ergy. So the dynamics of the fireball crucially depends on 
the baryon content of it. But irrespective of it, the bary- 
onic load has to be very small, otherwise, the expansion 
of the fireball will be Newtonian, which is inconsistent 
with the present observations. 

As discussed above, core collapse of massive stars and 
merger of binary compact objects (neutron star-neutron 
star, neutron star-black hole) as possible progenitors of 
the long and short GRBs respectively. In both the cases 
neutrons are the main ejected materials from the cen- 
tral engine. So it seems obvious that, the outflow from 
the central engine, which later on give rise to GRBs, 
has neutron as a principal component. The presence 
of neutron component in GRBs was first proposed by 
Derishev et al,[9| and later on it was showed that the 
presence of neutron component in the GRBs are prac- 
tically inevitable and its presence drastically changes 
the dynamics of the fireball. The possible role of neu- 
tron in observed emission has been discussed by many 
authors [l^, [ll|, [H, HI]. The dynamical decoupfing of 
neutron from the rest of the relativistic shell will give 
rise to inelastic n-p scattering and lead to emission of 
observable multi-GeV neutrinos [T^. ITsj . 

It is believed that, baryon number is not a good sym- 
metry of nature, due to the fact that, our universe is mat- 
ter dominant and this requires baryon number violating 
interactions. Thus the breaking of exact-baryon num- 
ber conservation in unified theories of the fundamental 
gauge interactions^]^ has lead to searches for proton de- 
cay and nn oscillation, but so far with no success 15 1. 
The nn oscillation occurs due to AB = 2 transition. 
The present experimental limit for unbound neutrons is 
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Tn-a ~ 10^ S0. 

The neutron oscillation in free space and its effect on 
cosmic rays has been discussed in Ref.^l?]. As a phe- 
nomenological implication of nn oscillation, it has been 
suggested that, the excess of sub-GeV anti-protons over 
protons in the galactic cosmic rays is due to the nn os- 
cillation in supernova explosions, which are site for large 
neutron excess[18j. Also it has been argued that, nn 
oscillation put a lower limit on the primordial magnetic 
field when applied to the neutron rich environment of the 
nucleosynthesis erap^. 

Here we propose that, if nn oscillations take place 
within the neutron rich GRB fireball, 19 — 38 GeV neu- 
trinos will be produced prior to the production of 5- 
10 GeV neutrinos due to dynamical neutron decoupling 
from the rest of the fireball as predicted by Bahcall and 
MeszarosfT^. By using earth as the detector and assum- 
ing the magnetic field in the jet outflow from the central 
engine to be '-^ 10^^ G, about 4 events can be recorded 
per year. The future international project Extreme Uni- 
verse Space Observatory (EUSO) may be able to detect 
these neutrinos [2Q]. 



II. nn OSCILLATION AND PHYSICS OF GRBS 

In the rest frame, neutron has a mean life time of Tp — 
888.6 ± 3.5 s. In the comoving frame, when the neutron 
has a Lorentz factor r„, the mean decay radius is rp = 
TnTp. So due to beta decay, at a distance r from the 
source, the neutron number is given by 

= A^„,oe-'^/'^^ (1) 

where Nnfi is original number of neutron within a ra- 
dius r ^ rp ^ 10^ cm. In the presence of a magnetic 
field, the neutron energy levels are split by an amount 
/S.E = gfiB and this is responsible for the oscillation 
of neutron to anti-neutron, where g — —1.91 is the 
anomalous magnetic moment of neutron, fi = eH/2mpC — 
3.152 X 10-12 gy Q-i and B is the magnetic field in 
Gauss. Due to nn oscillation, the number of anti-neutron 
is given bvp^ 

where Sm = Tnn^^ and here we take Tnn — lO^s. The 
number of n is inversely proportional to the square of 
the magnetic field. So for large magnetic field this is 
very much suppressed and for small one it is enhanced. 

In the fireball, a substantial fraction of the baryon ki- 
netic energy is transfered to a non thermal population 
of electrons through Fermi acceleration at the shock and 
these accelerated electrons cool via synchrotron emission 
and/or inverse Compton scattering to produce the ob- 
served prompt (due to internal shocks) and afterglow 
(due to external shocks) emission. For the synchrotron 



emission, strong magnetic field is needed to fit the ob- 
served data. But it is difficult to estimate the strength 
of the magnetic field from the first principle. One would 
expect large magnetic field if the progenitors are highly 
magnetized, for example, magnetars with B ~ 10^^ G. 
Also a relatively small pre-existing magnetic field can be 
amplified due to turbulent dynamo mechanism, compres- 
sion or shearing. Despite all these, there is no satisfactory 
explanation for the existence of strong magnetic field in 
the fireball. Also even if some magnetic flux is carried by 
the outflow, it will decrease due to the expansion of the 
flreball at a larger radius. Recently, it has been proposed 
that, the emission in the GRBs can also be explained 
through Compton-drag process and magnetic field is not 
necessary for the production of gamma rays [U [l^, 11^ . 
Because of the above uncertainties, here we consider the 
magnetic field as a parameter. 

By combining both the neutron decay and the nn oscil- 
lation, the number of anti-neutrons at a distance r from 
the source are given by 

where ^ is the neutron to proton ratio at the source and 
for simplicity we take ^ ~ 1 and the dimensionlcss en- 
tropy 7/ = E/M ^ 10^ — 10'^, where E is the initial ra- 
diation energy produced due to e*,7 and M is the rest 
mass of the fireball due to baryon load. 

The protons, neutrons and electrons are coupled to 
the thermal radiation in the expanding jet outflow from 
the central engine until the Compton scattering time 
scale t'rpf^ ~ (n'paTh)^^ and the elastic n-p scattering 
time scale ~ ("pCr„p)-i are shorter than the co- 
moving plasma expansion time scale t'^^p ^ r /T , where 
(Jnp ~ 3 X 10"^^ cm? is the n-p elastic scattering cross 
section, n'^ is the comoving number density of protons 
in the flreball. In fact when all the components in the 
jet are coupled, they have a common bulk Lorentz factor 
r ^ 300 Ref.[l2l- The neutrons and protons are coupled 
until the opacity T„p > 1, which corresponds to the r„p 
radius 
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(4) 

So the radius r for the n-p to be coupled should be smaller 
than r„p, i.e. r < r„p. Also the mean beta decay radius 
is rp ^ 10^^ cm >> r. So exp{—r/rp) ~ 1 and we can 
comfortably neglect the effect of beta decay in Eq.® 
when we are in the regime r < Vnp < rp. 

In the jet frame protons and neutrons have the same 
Lorentz factor F and the nn oscillation takes place with in 
the equilibrium plasma. The produced anti-neutrons can 
annihilate with the background neutrons to produce pi- 
ous and also can have elastic scattering (but not inelastic 
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scattering) with the background protons. The beta de- 
cay of anti-neutrons will be suppressed due to the longer 
life time of the anti-neutrons in the comoving jet plasma 
unless the decay takes place in the regime r < r„p < rp. 
We consider the following annihilation processes within 
the comoving plasma, 



n + n —>■ TT 
n + n —f TT^ 
n + n 



^0 



(5) 



The 7r° will decay to 2j and the charge pions will decay 



■-|-i/e(i/e)-|-i^^(i/^)-|-i/^(j/^). (6) 



In Eq.([5]) the energy carried by each pion in the cm. 
frame is 1.88 GeV/k.^ , where A:^ is the multiplicity of 
pion. The produced pions are relativistic and in the 
first step of Eq.(l6]), the charged pions will decay to 
/i* -t- v^{v^). In these decays, the average energy carried 
by muon is ~ 80% of the pion energy and rest is carried 
by t';i(i?/j)[23]- In the second step of the decay muon will 
decay as shown in Eq.Q and the average energy car- 
ried by each particle is about 1/3 of the muon energy. 
So on an average the energy carried by each neutrino or 
anti-ncutrino due to muon decay is ^ ~ Q.bGeV/kT^ in 
the comoving frame. In the observer frame this will be 
translated to 



GeV 1 

— 0.5 — I -T- ; r. 

K (1 + z) 



(7) 



By normalizing the burst at redshift z = 1, we obtain 

^v.v — '0 — i 300 



(8) 



For the pion multiplicity — i and 2, this gives e^^p — 
25 GeV and 37.5 GeV respectively. Similarly the energy 
of the muon neutrinos due to pion decay tt^ + 
is 



56 GeV/k^ T. 



300 



1 



(9) 



which is 19 GeV and 28 GeV respectively for k^^ ~ 3 and 
2. So the neutrinos and anti-neutrinos produced due to 
nfi annihilation are having energies in the range 19 GeV 
to 38 GeV. 

The neutral pions will decay to photons and have en- 
ergies 0.3 GeV (first of Eq.®) and 0.47 GeV (last of 
Eq.([5])) respectively in the cm. frame. In the observer 
frame, these will be boosted to 45 GeV and 71 GeV re- 
spectively. As we have discussed earlier, these photons 
are produced at a radius r < r„p and r^^ > 1, they will 
degrade their energy by producing pairs and thermal- 
ize within the fireball medium. But the neutrinos and 
anti-neutrinos which are produced due to nn annihila- 
tion will stream away from the fireball much before the 



decoupling of neutrons. So these neutrinos are different 
from the one that are produced due to neutron decou- 
pling in the later stage of the fireball evolutionp^. The 
total amount of energy released due to nn annihilation is 



e„ft = 2m„iVfi ~ 6 X lO^'^ ^^53'7ioo 



1+e 



erg, 
(10) 

and a major fraction of it will be emitted in the form of 
neutrinos and anti-neutrinos. 



III. POSSIBLE NUMBER OF NEUTRINO 
EVENTS 

By considering both the long and short GRBs rate 
within a Hubble radius of Rt ^ lO^Rb5^^, the number 
of events per year in a detector with Nt target protons is 



where 



D = 2.64 X 10 



28 



^70^(1 ^ l/Vl + z) cm, 



(11) 



(12) 



is the proper distance out to red shift z with /lyo = 
Ho/70 km/s/Mpc and a^p - 5 x 10""*"r(l + cra^ 
is the total detection cross section per neutron. If we 
consider the whole earth as the detector with Nt ~ 10^^ 
and taking a small magnetic field B ^ 10^^ G in the jet 
outflow (in Eq.(I2])), the event rate is 



Ruu 



4 h'^QRb5 — — -/VfsiFsoo 
f?ioo 




2 + z 




year 



''•(13) 



So the expected rate of multi-GeV neutrino event is 
about 4 per year and the increase in the initial neutron 
to proton ratio ^ will slightly increase the event rate[26j. 
May be in future, the EUSO will be able to detect these 
neutrinos by emission of fluorescence light of nitrogen 
due to the extensive air showers through the interaction 
of the multi-GeV neutrinos with the atmosphere (20j. 



IV. CONCLUSIONS 

The production of neutrinos and anti-neutrinos due to 
nn annihilation is inversely proportional to the square 
of the magnetic field. If the magnetic field in the fire- 
ball is very weak, the number of anti-neutrons will be 
more as can be seen from Eq.(I2|). The produced anti- 
neutrons will annhilate with the background neutrons 
through the processes in Eq.(l5]) to produce neutrinos, 
anti-neutrons and photons. A large fraction of the an- 
nihilation energy will stream away in the form of neu- 
trinos, thus reducing the energy content of the fireball. 
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For example if the magnetic field in the fireball is very 
small {B ~ 10~^^ G), the amount of energy release due 
to nh annihilation is e„fi ~ 10^^ erg. Apart from the 
nn annihilation, if beta decay of anti-neutrons take place 
within the region r < Vnp with a very small magnetic 
field, the excess amount of baryon, anti-baryon annihila- 
tion will result in decreasing the baryon content of the 
fireball further and making it cleaner, which can prob- 
ably explain the almost baryon free environment of the 
fireball. On the other hand, strong magnetic field will 
suppress the nn oscillation and ultimately produce very 
less number of neutrinos and anti-neutrinos. So if the nn 
oscillation exist in GRB fireballs, magnetic field can not 
be arbitrarily weak or strong. The produced neutrinos 
and anti-neutrinos are unique signature of nh oscillation 
as well as the presence of magnetic field in the GRB fire- 
balls. It is unique in the sense that, there is no other 
process which can produce these neutrinos in between 
the production point (the source at r = rp ) and r — r^p- 
Flux of these neutrinos depends on the original neutron 



contain as well as the magnetic field in the fireball. Due 
to very low flux, these neutrinos will be very difficult to 
detect. The future project EUSO may be able to detect 
these neutrinos. However, if at all detected, it will be 
before the 5-10 GeV neutrinos due to inelastic scattering 
of decoupling neutrons with the protons in the fireball. 
The observed time difference between the nh neutrinos 
and the 5-10 GeV neutrinos is At ~ Vnp/T^ ~ 4 x 10^^ s. 
Although this time difference is very small, the energies of 
the neutrinos are very different. So the neutrinos within 
the energy range 19-38 GeV have different physics (pro- 
duction mechanism) as compared to 5-10 GeV neutrinos. 
Observation of these multi-GeV neutrinos will not only 
tell more about the nature of the GRB progenitors but 
also shed more light on the physics beyond the standard 
model. 
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